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In order to probe with high precision the tails of the ground-state energy distribution of disordered 
spin systems, Korner, Katzgraber and Hartmann have recently proposed an importance-sampling 
Monte-Carlo Markov chain in the disorder. In this paper, we combine their Monte-Carlo procedure 
in the disorder with exact transfer matrix calculations in each sample to measure the negative 
tail of ground state energy distribution Pd(Eo) for the directed polymer in a random medium of 
dimension d = 1,2, 3. In d = 1, we check the validity of the algorithm by a direct comparison 
with the exact result, namely the Tracy- Widom distribution. In dimensions d = 2 and d — 3, we 
measure the negative tail up to ten standard deviations, which correspond to probabilities of order 
Pd(Eo) ~ 10~ 22 . Our results are in agreement with Zhang's argument, stating that the negative tail 
exponent n(d) of the asymptotic behavior InPd(Eo) ~ — \Eo\ v ^ as Eq — > — oo is directly related to 
the fluctuation exponent 9(d) ( which governs the fluctuations AEq(L) ~ L e ™ of the ground state 
energy Eo for polymers of length L) via the simple formula T](d) = 1/(1 — 9(d)). Along the paper, 
we comment on the similarities and differences with spin-glasses. 



I. INTRODUCTION 



Since the ground-state energy Eq of a disordered sample is the minimal energy among the energies of all possible 
configurations, the study of its distribution belongs to the field of extreme value statistics. Whereas the case of 
independent random variables is well classified in three universality classes 0] , the problem for the correlated energies 
within a disordered sample remains open and has been the subject of many recent studies. The interest lies both 

(i) in the scaling behavior of the average Eq v (L) and the standard deviation AEq(L) with the size L 

(ii) in the asymptotic distribution P(x) of the rescaled variable x — (Eq — Eq v (£)) / AEq(L) in the limit L — > oo 

In this introduction, we first recall what is known in the field of spin-glasses, before focusing on the directed polymer 
model. 



A. Ground state energy distribution in spin-glasses 

For spin-glasses in dimension d, let us consider samples containing N — L d , where where L denotes the linear size, 
and follow the notations of Ref. 0]. The ' shift exponent' S governs the correction to extensivity of the averaged 
value 

Eq V (L) ~ L d e a + L ° ei + ... = Ne + N e °' d e 1 + ... (2) 

Within the droplet theory 4, 5], this shift exponent S coincides with the domain wall exponent Odw and with the 
droplet exponent 6 of low energy excitations. The ' fluctuation exponent' Of governs the growth of the standard 
deviation 

AE (L) - L fe 2 = N e f' d e 2 (3) 

In any finite dimension d, it has been proven that the fluctuation exponent is Of = d/2 raL Accordingly, the rescaled 
distribution P(x) of Eq. Q was numerically found to be Gaussian in d = 2 and d = 3 suggesting some Central 
Limit theorem. On the contrary, in mean-field spin-glasses, the width does not grows as N 1 ^ 2 and the distribution 
is not Gaussian. In the Random Energy Model 0, the width remains finite AEq(N) ~ 0(1) and the distribution is 
the Gumbel distribution In the Sherrington-Kirpatrick model, the width grows as AEq(N) ~ N 1 / 4 according to 
some theoretical arguments 0, and numerics 0, ll(| , and the distribution is clearly asymmetric 0, ll(| . Finally for 
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the one dimensional disordered spin chain with power-law interactions that allows to interpolate between effectively 
finite-dimensional and mean-field models, the transition between short-range and infinite-range behaviors corresponds 
to the Gaussian-non Gaussian transition for the ground state energy . 



B. Ground state energy for the directed polymer 



The directed polymer model in 1 + d dimensions, is defined by the following recursion for the partition function 

2d 



(4) 



The bond-energies e^f + ej,r) are random independent variables, drawn with the Gaussian distribution 



(5) 



This model has attracted a lot of attention for two main reasons : (i) it is directly related to non-equilibrium 
jroperties of g rowth models |12| (ii) as a disordered system, it presents some similarities with the spin-glass physics 
At low temperature, there exists a disorder dominated phase, where the order parameter is an 




The probability distribution of the ground state energy Eq is expected to follows the scaling form of Eq. (|TJ . In 
contrast with spin-glasses where the shift exponent S (Eq. 0) and the fluctuation exponent Of (Eq. 0) are different, 
there is a single exponent 0(d) that governs both the correction to extensivity of the average Eq v (L) and the width 
AE (L) 



E$ V (L) 
AE (L) 



Le + L eid) e 1 + ... 
L e ^e 2 + ... 



(6) 
(7) 



This exponent also governs the statistics of low excitations within the droplet theory p|, as confirmed numerically 
|l6j . This exponent is exactly known in one-dimension [rH fl8l IT^. l20j 

0(d=l) = l/3 (8) 

and has been numerically measured in dimensions d — 2,3,4,5 0, |2^, |2^, 0, HH 

9(d = 2) - 0.244 (9) 
9(d = 3) -0.186 (10) 



For the mean-field version on the Cayley tree, the exponent vanishes 0(d = oo) = 13, 26], with a width of order 
0(1) for the probability distribution, but with a non random O(lnL) correction to the extensive term e$L in the 
averaged value [l3j . 

The rescaled distribution Pd is exactly known in d = 1 and is related to Tracy- Widom distributions of the largest 
eigenvalue of random matrices ensembles ^3H3;0|- O n the Cayley tree, the rescaled distribution was found to be 
non universal and to depend on the disorder distribution |26|. 



C. Numerical measure of the ground state energy distribution 

The numerical measure of the ground state energy distribution is usually done by a simple sampling procedure, 
where the histogram of the energies of independent samples are collected. However recently, Korner Katzgraber and 
Hartmann have proposed an importance-sampling Monte-Carlo algorithm in the disorder, which allows to measure 
much more precisely the tails of the distribution. In the case of the Sherrington-Kirkpatrick model of spin-glasses, 
this procedure was used to measure the negative tail on systems of size N < 128 Q] up to x > — 15 corresponding 
to probabilities P(x) > 10~ 18 (see Eq. QJ, whereas the simple sampling procedure cannot go beyond x > — 5 
corresponding to P(x) > 10~ 4 |10| . 

For the directed polymer in dimensions d — 2 and d = 3, the rescaled distribution Pd has been numerically measured 
via simple sampling i n|27| with results in the region x > —5. In this paper, we use the importance sampling algorithm 
recently proposed in [l| to measure precisely the negative tail of the probability distribution up to x > —10. 
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The paper is organized as follows. In Section |HJ we recall Zhang's argument 0] that relates the decay of the 
rescaled distribution P d to the fluctuation exponent Of. In Section UTTl we describe the Monte-Carlo procedure in the 
disorder proposed in [jj and mention the specific choices for the application to the directed polymer model. In Section 
IIVI we show the validity of the procedure in d = 1 via the direct comparison with the exactly known distribution 
(Tracy- Widom). Finally in Sections Ivl and Ivl! we present our results for d = 2 and d — 3 respectively. We present 
our conclusions in Section fVIII 



II. ZHANG'S ARGUMENT FOR THE NEGATIVE TAIL EXPONENT 



A. Distribution of the free-energy in the low temperature phase 

According to the droplet theory Q, the whole low temperature phase < T < T c is governed by a zero-temperature 
fixed point. In particular, at T < T c , the droplet exponent 9(d) governs the width AF(L,T) and the correction to 
extensivity of the average F av (L, T) 

AF(L, T) ~ L e ^f 2 (T) + ... (11) 
F av (L,T) ~Lf (T) + L e Wf 1 (T) + ... (12) 

and the rescaled probability distribution of the free-energy coincides with the rescaled distribution Pd describing the 
ground-state energy distribution Q 

W L, T) * P d (, = F l F p l {L T P ) (13) 

as recently checked numerically using simple sampling [28l | . 



B. Zhang's argument for the directed polymer 

In finite dimensions d > 1, the rescaled distribution P c i is not known but there exists a simple argument due to 
Zhang 0| that allows to determine the exponent r\ of the negative tail of the free energy distribution 

P d (x -> -oo) ~ e - c ^ v<d) (14) 

If r/(d) > 0, the moments of the partition function can be evaluated by the saddle-point method, with a saddle value 
F* lying in the negative tail l|14|) 



Zf = J dFV L (F, L)e~ nf)F ~ / dFe- c \^) e~ n ? F ~ e Hn)L^r (15) 

Since for positive integer n, these moments of the partition function can be formulated in terms of the iteration of 
some transfer matrix, they have to diverge exponentially in L with some Lyapunov exponent. As a consequence, the 
exponent 77(d) of the negative tail (|14fl is not a free parameter, but is fixed by the value of the fluctuation exponent 

^ - i-W (i6) 

In dimension d = 1 where the droplet exponent is exactly known 6(d = 1) = 1/3 (Eq. |SJ, this yields the negative tail 
exponent 

V(d=l) = l (17) 

in agreement with the exact Tracy- Widom distributions [T^. I20ll2l| . In dimensions d = 2 and d = 3, the numerical 
estimates of the droplet exponents (Eq. I10|l yield the following predictions 

r)(d = 2) -1.32 

rj(d = 3) -1.23 (18) 

These predictions have been tested numerically in |27j using simple sampling that do not allow to have data far in 
the tails. In the following, we will use the importance sampling Monte-Carlo method in the disorder to probe the 
negative tail more precisely. 
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C. Zhang's argument for spin-glasses 

To the best of our knowledge, Zhang's argument seems to be applied only in the context of directed polymers |12| . 
whereas it can be applied for other kinds of disordered systems since it is only based on scaling argument within a 
saddle-point approximation in the large L limit (Eq. I15|) . It is thus interesting to describe now its implications in the 
field of spin-glasses. 

For spin- models in finite dimensions, the fluctuations of free energies over the samples scale instead as 

\^F L ] samples ~ L d ^ (19) 

at any temperature as proven in jy]. This scaling simply reflects the Central-Limit fluctuations of the L d disorder 
variables defining the sample. ( The directed polymer escapes from these normal fluctuations because it is a one- 
dimensional path living inal + d disordered sample : each configuration of the polymer only sees L random variables 
among the L 1+d disorder variables that define the sample, and the polymer can 'choose' the random variables it sees.) 

Repeating Zhang argument in this case (|19|l yields for the negative tail exponent rj(d) = 2. This is in agreement with 
the recent numerical studies 0, ITU ] that find a Gaussian distribution in finite dimensions and in the one dimensional 
Ising spin-glass with long range interactions in the non-mean-field regime. 

On the contrary, for the Sherrington-Kirkpatrick model 0, |jj ll(| , the probability distribution of the ground state 
is found to be asymmetric, and has been fitted with generalized Gumbel distribution [ltllOf. 

However, if one repeats Zhang argument for the SK model with the measured fluctuation exponent Of ~ 0.235 [Tel ] 
for the width AEq (N) ~ N e f , one obtains the negative tail exponent 

nsK = y^-jj ~ i.3 (20) 

If the value of the fluctuation exponent is exactly Of — 1/4 as suggested by some theoretical arguments 

dEl, the 

negative tail exponent would be tjsk = 4/3. 

This could explain why the fit with generalized Gumbel distributions whose negative tail is a simple exponential 
g -m\x\ w jt n exponent rj = 1 and coefficient m leads to increasing effective values of m when the range over which 
the tail is measured grows : the fit with simple scaling data oni> —6 leads to m ~ 6 |10| . whereas the importance 
scaling data on x > —15 leads to a completely different estimate m ~ 11 0. 



III. DESCRIPTION OF THE IMPORTANCE-SAMPLING MONTE-CARLO ALGORITHM IN THE 

DISORDER 

In Ref. |H, a procedure based on an importance-sampling Monte-Carlo algorithm in the disorder was proposed to 
probe with high precision the tails of the ground-state energy distribution of disordered systems, and was applied 
for the Sherrington-Kirpatrick mean-field Ising spin-glass, where probabilities up to 10~ 18 could be measured. In 
this Section, we summarize their method which can be divided in three steps. For each step we mention the specific 
choices we have made to apply it to the directed polymer model. 



A. Simple sampling 

A disorder configuration will be denoted by T>, and its ground state energy E(T>). For the directed polymer, the 
ground state energy can be computed via transfer matrix. A simple sampling numerical estimation P S impie {E) of the 
ground state energy distribution P(E) consists in drawing n s independent disordered samples T>\, ...T> n ^, in computing 
the corresponding ground state energies E(T>i), E{T> n3 ), and in constructing the histogram 



P simp i e {E) = —Y,b{E-E{V i )) (21) 

This histogram is very useful to measure the distribution P(E) where P{E) ^> — , but gives no information on the 

tails where P{E) < since no events are found. 

As an example, we show on Fig. ^ l ne results we have obtained recently via simple sampling for d = 1,2,3 
respectively [jjj : whereas the core of the distribution is well measured, the tails suffer from statistic fluctuations 



FIG. 1: (Color online) Logarithmic plot of the rescaled probability distribution Pd(x) (Eq. as measured via simple 

sampling for (a) L = 12, 50, 200, 800 in d = 1 ; (b) L = 10, 40, 160 in d = 2; (c) L = 6, 18, 36 in d = 3 



as soon as the probability becomes too small. Moreover, if one chooses to make the same CPU effort on all sizes, 
the number of samples rapidly decay with the size L, so that the data for the tails are less and less precise as L 
grows. Since one is interested into the asymptotic regime L — > oo, the correct measure of the tails quickly becomes 
intractable within the simple sampling procedure. 

This is why a correct measure of the tails requires the use of some importance sampling, as stressed in However, 
the simple sampling study is the first necessary step within the present method, for three reasons : 

(i) the simple sampling results are needed to construct the guiding function of the importance sampling measure^ 
as described below 

(ii) the simple sampling results give accurate results for the average value Eq v (L) and the standard deviation AEq(L), 
that do not have to be measured via importance sampling 0. In particular, this allows to work on a finite box 
[xmin, x-max] for the rescaled variable and to choose freely the boundaries of the box, for instance x max = — 1. to 
concentrate on the negative tail, as will be done below for the directed polymer. 

(iii) finally, the simple sampling results allow to check the validity of the importance sampling measures on the core 
of the distribution where the simple sampling results are sufficiently precise. 

B. Construction of a guiding function G(E) from the simple sampling result P S im P ie(E) 




FIG. 2: (Color online) Fits of the GOE Tracy- Widom distribution In P^w ( x ) by generalized Gumbel distributions lng m (x) 
(Eq. 1221 on various intervals : (a) fit on [-5,3.7] corresponding to InP > —10 with m = 12.93 ; (b) fit on [-8,5] corresponding 
to InP > -20 with m = 14.71 (c) fit on [-10,6] corresponding to InP > -30 with m = 15.92 ; 

The simple sampling result P S impie{E) exists in the range in E where P s impie(E) > l/n s , whereas the guiding 
function G(E) needed for the importance sampling below has to be defined in the tails where P S im P ie(E) < l/n s . 
The guiding function G(E) should be in some sense the 'best' extrapolation of the data P simp i e (E). The proposal of 
is to define G(E) as the best fit of P S impie(E) within the one parameter family of generalized Gumbel distribution 
g m (x), which reads for the normalization conditions < x >— and < x 2 >= 1 
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where (3{m) — , 1 and aim) = — /3(m) ( p/"V — In m J . The usual Gumbel distribution corresponds 

,/ r"tm) f r'( TO ) \ 2 \ (m) / 

to to = 1, whereas the Gaussian can be formally recovered in the limit to — > oo. This choice was motivated by 
the numerical finding that the rescaled probability distribution as measured via simple sampling could be fitted 
extremely well by a generalized Gumbel distribution g m (x) of parameter m ~ 6 0, . It turns out that many recent 
studies in various contexts have found that asymmetric distributions could be extremely well fitted by generalized 
Gumbel distributions with various non-integer values of to [2^, [3(3, [^J . This has motivated theoretical studies to 
understand the origin of this type of distribution 33] . However as discussed in ref. |29j , it is empirically known that 
probability distribution functions (PDF) with the same first four moments approximately coincide over the range of 
a few standard deviation which is precisely the range of numerical -or experimental- data. So generalized Gumbel 
PDF's, with arbitrary to, should not be considered more than a convenient one parameter fit. To demonstrate clearly 
how misleading these fits can be, we show on Fig [21 how the Tracy- Widom GOE distribution P^w E (which represents 
the exact rescaled distribution for the directed polymer model in d = 1 0, [2(1 EI) can be fitted by generalized 
Gumbel distributions on the three ranges P > 1CP 10 , P > 1CP 20 and P > 10 -30 : the best fit corresponds to 
increasing values of the parameter to. For the first range P > 10~ 10 , the found fit is 'perfect', whereas a slight 
difference begin to appear in the negative tail as the range grows. Moreover, the Tracy- Widom distribution is known 
to have the following asymptotic behavior 

P^ E (x) ~ e-^l" 1 with 77! = % (23) 

x — > — oo Z 

whereas the generalized Gumbel distribution have for for any m an exponential tail with exponent ry = 1 and coefficient 
m 

g m (x) ~ e~ m W (24) 

x — > — oo 

This explains why the effective m of the best fit grows with the range. In conclusion, whenever the fit of the core of 
the distribution leads to an effective to which grows with the range, as in the SK model where to ~ 6 and to ~ 11 
were found depending on the range 0,113, the PDF is probably not a generalized Gumbel distribution, but is likely 
to have a negative tail exponent r\ > 1 (as already suggested around Eq. 1201 using Zhang's argument ). And if one 
focuses on the negative tail, it is clear that the fit with a simple exponential l|24|l is very restrictive. 

As a consequence, in the following where we focus on the negative tail x < — 1 for the directed polymer, we have 
chosen not to work with generalized Gumbel distribution, but to construct a guiding function G(E) which fit the 
simple sampling data and whose leading behavior involves the negative tail exponent 77^ as obtained from Zhang 
argument ( see Eas ll6l |17I and lT%|) . In practice, we have found convenient to work in d = 2 and d — 3 on the range 
x e [—10, —1] with some guiding function Gd{x) of the form 

\nG d {x) = ao- a\\x\ Vd + a,2\n\x\ (25) 

where the three parameters a,i{d) were chosen to fit best the simple sampling data. 



C. Importance sampling with the guiding function G(E) 

The importance sampling Monte Carlo algorithm proposed in y| is defined by the following Markov chain : 

(1) From the current disorder configuration T>i, construct a candidate V for the next disorder configuration T>i + \ 
by replacing a subset of T>i chosen at random with new values drawn with the original disorder distribution. For a 
spin model of N spins, this subset can be for instance a single bond chosen at random, or all bonds connected to 
a site chosen at random, so that the proposed change in the ground state energy is of order O(l) with respect to a 
value of order Eq(N) = Neo + i.e. its relative order of magnitude is of order 1/N [J|. For the directed polymer 
studied here, we have chosen for this subset the energies of a whole time-slice, i.e. all the disorder variables seen by 
a given monomer. Then the proposed change in the ground state energy is of order 0(1) with respect to a value of 
order Eo(N) = Neo + •■■ as in spin models. 

(2) Calculate the new ground state energy E(V') and compare it with the previous ground state energy E(T>i) 
using the guiding function G(E) : set T>i + \ = V with probability 



V accept (T>'\T>i 



G{E{V t )) 
G{E{V')Y 



(26) 



and set T>i + \ — T>i otherwise. 
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This Markov chain is expected to converge towards a stationary state where a disorder configuration V is visited 
with probability cx \/G(E(T>)). The stationary probability to visit a disorder configuration with energy E is now 
given by the ratio 

Rstationary(E) = QTj^j ffiO 

If the guiding function G(E) were the exact P(E), this would corresponds to a flat-histogram sampling of P{E). If 
G{E) is just a reasonable extrapolation of the simple sampling result P S i mp i e (E), one expect to measure nevertheless 
much better the tails of P(E). 

(3) Measurements from the Monte-Carlo procedure : since successive configurations visited by a Monte-Carlo 
algorithm are not independent, one usually keeps only decorrelated configurations for the numerical measure 
Rimportance(E) of the theoretical stationary solution R s tationary(E). This means in practice that one should first 
estimate some typical correlation time r and use only every rth configuration 

j mi 

RtlportancJE) = — £ 6 (E - E(V l+JT )) (28) 
1 3=1 

where the number mj of measured is simply the ratio mj = — of the total number T of Monte-Carlo iterations by 
the correlation time r. For instance in [l|, the time r was chosen to be r — Ar e where r e is the time where the 
autocorrelation of the ground state energy 

nl .\ _ < EiEi+t > - < Ei >< E i+t > . , 

C[t) ~ <E?>-< % >s (29) 

decays to 1/e. For the SK model with 16 < N < 128 spins, the autocorrelation time was found to be of order of 
400-700 MC steps Q. 

For the directed polymer, we actually find that the histograms Rim por tance(E) obtained for r = 1 and t r e 
coincide, except that the histograms with large r contain more noise since they are built out of less events. From a 
theoretical point of view, one can justify this finding as follows : if the total Monte-Carlo time T is much bigger than 
the typical time t cross to cross the interval [E m i n ,E max ], then the average with respect to the stationary measure 
should be equivalent to the time average of the Monte-Carlo procedure where all times are kept 

r i T 

/ dE f i^E^j ^stationary 

(E) = -J2f( E (t)) fOT T^>t cross (30) 

Indeed for a free random walk in a finite box, it seems clear that one obtains the flat histogram via measuring the 
positions at all times, instead of throwing away most of the times to have independence between two consecutive 
measures. The quality of the convergence towards the stationary distribution then depends on the number 

T 

Across ~ (^-Q 



of crossings of the interval [E m i n , E max ] during the total number T of the Monte-Carlo, which should be large enough 

Across 1- 



D. Summary of the procedure used for the directed polymer 

In the following sections, we will present the results for the ground state energy distribution obtained by combining 

(i) the Monte-Carlo procedure in the disorder discussed above 

(ii) the transfer matrix calculation of the ground state energy in each sample with a free boundary condition for the 
end polymer. 

We have chosen to focus on the negative tail, by working on the same finite box x € [x m i n ,x max ] in terms of the 
rescaled variable x (Eq. ^| for all sizes L. We now present our results for d — 1,2,3 respectively. 
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IV. MEASURE OF THE GROUND STATE ENERGY DISTRIBUTION IN d = 1 



In d = 1, the exact rescaled distribution of the ground state energy is exactly known and corresponds to the 
Tracy- Widom GOE distribution P<f° E [H ifl |2l| if the last monomer is free, the case we consider here. (It would 
be the Tracy- Widom GUE distribution P^w E ^ the ^ as * monomer were fixed at the origin) . We use this exact result 
to check the validity of the Monte-Carlo procedure in the disorder and to describe its main properties. 



A. Numerical details 



In dimension d = 1, we have chosen to work on the interval [x m i n , x ma x\ = [ — H-,— 1.] for the rescaled variable 
x (Eq. nj, i.e. to probe the negative tail up to probabilities of order Pi(x) > 10~ 32 . We now give the sizes L we 
have studied, together with the standart deviation AEq(L) measured by simple sampling and used in the rescaling 
of Eq. ^ (the averaged values Eq v (L) can be found in our previous work 16] ), the corresponding number Tl of 
Monte-Carlo iterations, the acceptation rate T acc (L) of Monte-Carlo moves, and the number n cross (L) of crossings of 
the box [x min ,x max ] = [-11., -1.]. 

L = 50,100,200,400,800,1600 (32) 
AE (L) - 3.12,3.98,5.04,6.36,8.04,10.11 (33) 
T L = 33.10 8 ,85.10 7 ,225.10 6 ,57.10 6 ,7.10 6 ,4.10 6 (34) 
Tacc(L) - 0.54,0.64,0.72,0.76,0.82,0.86 (35) 
n cr oss{L) - 98.10 4 ,224.10 3 ,45.10 3 ,8.10 3 ,670,230 (36) 

As L grows, the proposed Monte-Carlo moves Ax in the rescaled variable x (Eq. ^) are smaller : this is why both 
the acceptation rate T acc (L) and the crossing time t cr oss{L) (Eq. 131 [I also grows with L. The final result is that the 
number of crossing n cross (L) decays with L, and since it should remain large enough to obtain a good measure (Eq. 
1310 . this number fixes the maximal size that can be correctly studied. 



B. Properties of the Monte-Carlo process in the disorder 



3e+08 



1e+08 




1e+05 



50000 - 




FIG. 3: (Color online) Monte-Carlo procedure to measure the negative tail on x G [—11,-1] in d = 1 : histograms of the 
proposed and accepted Monte-Carlo changes Ax in the disorder (a) for L = 25 where the acceptation rate is r aC c ~ 0.425 (b) 
for L — 1600 where the acceptation rate is r acc ~ 0.86 . 



We show on Fig. [31 the histograms of the proposed and accepted Monte-Carlo changes in the disorder, for L = 25 
and for L = 1600 respectively. The proposed changes are biased towards Ax > 0, because here, in the negative tail, a 
Monte-Carlo step Ax > corresponds to a move where the probability P(x) is bigger. The histogram of the accepted 
moves is on the contrary almost symmetric around Ax — in order to generate a non-biased random walk. For 
Ax < 0, the two histograms almost coincide, i.e. a move Aa; < is almost always accepted. As the size L grows, the 
proposed moves in the relative variable x are smaller, and as a consequence, the acceptation rate grows with L. 
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t t 

FIG. 4: (Color online) Monte-Carlo procedure to measure the negative tail on x 6 [—11, —1] in d = 1 : process x(t) during the 
first 1 < t < 10000 Monte-Carlo iterations (a) for L = 50 (b) for L = 200 



The resulting process x(t) are shown on Fig. Q]for the first 1 < t < 10000 Monte Carlo iterations, for L = 50 and 
L = 200 respectively. The time t cross needed to cross the interval [x m i n ,x max ] = [—11, —1] grows with L. 



C. Convergence towards the exact Tracy- Widom distribution 



0.15 



0.13 - 



0.09 



0.07 



0.05 




(a) 



In P (x) 
1 



-10 - 



-20 



-30 - 



-40 





X \ - 

€ V 


(b) 





-15 



-10 



FIG. 5: (Color online) Monte-Carlo procedure to measure the negative tail on x £ [—11, —1] in d = 1 : (a) relative histogram 
Pl(x) — Pl(x) I Ptw E \ x ) with respect to the exact guiding function : convergence towards the flat histogram as L grows 
: L = 200 (dashed line), 400,800, 1600 (thick line), (b) logarithmic plot of the negative tail of the probability distribution 
Pi(x), as compared to simple sampling result for L = 1600. The exact Tracy-Widom distribution is also shown (thin line) to 
demonstrate the validity of the Monte-Carlo procedure. 



On Fig. [U(a), we show the relative histogram Pl(x) / P^r E (x) of the measured Pl{x) via the Monte Carlo procedure 
as L grows with respect to the Tracy-Widom GOE distribution that represents the asymptotic exact result for L — > oo 
: these relative histograms becomes natter as L grows. 

On Fig. 0(b), we show for comparison : 

(i) the simple sampling histogram for L = 1600 

(ii) the importance sampling measure of the tail on x G [—11,-1] for L = 1600 

(iii) the exact Tracy-Widom GOE distribution. 

Our conclusion is thus that the Monte-Carlo in the disorder is a very efficient method to probe accurately the tails, 
since they allow to reproduce the exact result on the range x e [—11, —1] for sizes up to L = 1600. 
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D. Extraction of the negative tail exponent value 



Let us now make some comments on the extraction of the negative tail exponent value. The exact Tracy- Widom 
GOE distribution P^^ E {x) has for negative exponent rji = 3/2. However, the following fits of this distribution 
P!f° E ( x ) on the finite range x G [—11,-1] give slightly larger values : 

(i) the fit of (In P^^ E {x)\ by a — b(— x) Vl containing three parameters yields r/i ~ 1.58 

(ii) the fit of (In Pj!{y E (x)) by a — b(—x) Vl + cln(— x) containing four parameters yields 771 ~ 1.54. 

This shows that the extracted value of the negative tail exponent from data on the finite range x G [—11,-1] is not 
very precise if there is not information on the subleading terms. 

Similarly in higher d below, we expect that the Monte-Carlo procedure gives very accurate data on the range where 
the tail is measured, but that the extraction of the negative tail exponent value suffers from some error directly related 
to the range that is probed. 



V. RESULTS FOR THE GROUND STATE ENERGY DISTRIBUTION IN d = 2 



A. Numerical details 



In dimension d = 2, we have chosen to work on the interval [x m i n ,x max ] = [—10., —1.] for the rescaled variable x 
(Eq. [IJ, i.e. to probe the negative tail up to probabilities of order P\{x) > 10 -23 . We now give the sizes L we have 
studied, together with the standart deviation AEq(L) measured by simple sampling and used in the rescaling of Eq. 
U(the averaged values E§ V (L) can be found in our previous work the corresponding number Tl of Monte-Carlo 
iterations, the acceptation rate T acc {L) of Monte-Carlo moves 

L = 20,40,80,120,160 (37) 

AE a (L) - 1.58,1.85,2.18,2.40,2.54 (38) 

T L = 125.10 7 ,27.10 7 ,47.10 6 ,34.10 5 ,64.10 4 ,21.10 4 (39) 

Tacc(L) - 0.24,0.28,0.38,0.47,0.5,0.54 (40) 



B. Monte-Carlo results 




2000 4000 6000 8000 10000 

t 



FIG. 6: (Color online) Monte-Carlo procedure to measure the negative tail on x G [—10, — 1] in d — 2 for L — 120 : (a) 
process x(t) during the first 1 < t < 10000 Monte-Carlo iterations (b) logarithmic plot of the negative tail of the probability 
distribution ^(a;), as compared to simple sampling result. 



On Fig. El (a), we show the process x(t) during the first 10 000 Monte-Carlo iterations for L = 120. On Fig. 03(b), 
we compare the importance sampling measure of the negative tail with respect to the simple sampling evaluation. 
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C. Negative tail exponent 7^=2 

From the point of view of the convergence in L towards a fixed distribution, we find that the negative tail measured 
for the two bigger sizes L = 120 and L = 160 nearly coincide on the whole interval [x m i n , x max ] = [—10., — 1.] under 
study ( whereas our results for the smaller sizes do not). 

As explained previously in Section TlV Dl for the case d = 1, the error on the estimated value of the negative tail 
exponent is due to the range [x m i n , x max ] = [— 10.,— 1.] over which the fits are made. As in Section llV Dl we have 
tried to fit our result for (lnJ^aO) as measured for the sizes L = 120 and L = 160 by the two following fits, with or 
without power-law corrections with respect to the leading exponential term : 

(i) the first fit a — b(— x) n2 containing three parameters yields 772 ~ 1.4 

(ii) the second fit by a — b(—x) m + cln(— x) containing four parameters yields 772 — 1.3. 

Our conclusion is thus that the extracted value of the negative tail exponent from our data on the finite range 
x G [—10, —1] is not very precise in the absence of information on the subleading terms, but is compatible with the 
value — 1.32 predicted by Zhang's argument (see Eqs. 1161 and IT8*|) . 

VI. RESULTS FOR THE GROUND STATE ENERGY DISTRIBUTION IN d = 3 

A. Numerical details 

In dimension d = 3, we have chosen to work on the interval [x m i n , x max ] = [— 10.,— 1.] for the rescaled variable 
x (Eq. nj, i.e. to probe the negative tail up to probabilities of order P\{x) > 10~ 21 . We now give the sizes L we 
have studied, together with the standart deviation AE (L) measured by simple sampling and used in the rescaling 
of Eq. ^ (the averaged values E§ V (L) can be found in our previous work 16] ), the corresponding number Tl of 
Monte-Carlo iterations, the acceptation rate r acc (L) of Monte-Carlo moves, and the number n cross (L) of crossings of 



the box [x min ,x max ] = [-10., -1.]. 

L = 12,24,36,48,60,72 (41) 

AE {L) - 1.15,1.30,1.39,1.46,1.52,1.55 (42) 

T L = 64.10 6 ,43.10 5 ,75.10 4 ,95.10 4 ,34.10 4 ,182.10 3 (43) 

TacciV) - 0.24,0.27,0.32,0.34,0.36,0.37 (44) 

n cross {L) - 20000,2400,300,522,166,94 (45) 



B. Monte-Carlo results 




FIG. 7: (Color online) Monte-Carlo procedure to measure the negative tail on x G [—10, —1] in d — 3 for L = 72 : (a) process 
x(t) during the first 1 < t < 10000 Monte-Carlo iterations (b) logarithmic plot of the negative tail of the probability distribution 
Ps(x), as compared to simple sampling result. 
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On Fig. 0(a), we show the process x(t) during the first 10 000 Monte-Carlo iterations for L = 72. On Fig. 0(b), 
we compare the importance sampling measure of the negative tail with respect to the simple sampling evaluation. 

C. Negative tail exponent rjd=s, 

As in Section HVDI we have tried to fit our result for (h\P^{x)) by the two following fits, with or without power-law 
corrections with respect to the leading exponential term : 

(i) the first fit a — b(—x) V3 containing three parameters yields 773 ~ 1.25 

(ii) the second fit by a — b(—x) rn + cln(— x) containing four parameters yields 773 ~ 1.15. 

Our conclusion is that the extracted value of the negative tail exponent from our data on the finite range x £ 
[—10, —1] is not very precise in the absence of information on the subleading terms, but is compatible with the value 
t]% = 1.23 predicted by Zhang's argument (see Eqs. [TBI and IT5jl . 

VII. CONCLUSION 

In this paper, we have adapted the importance-sampling method in the disorder proposed in |lj for spin-glasses, to 
measure with high precision the negative tail of the ground-state energy distribution P,i(E ) for the directed polymer 
in a random medium of dimension d = 1,2, 3. In e? = 1, we have checked the validity of the procedure by a direct 
comparison with the exact result, namely the Tracy- Widom GOE distribution. In dimensions d = 2 and d — 3, we 
have measured the negative tail up to P ~ 10~ 22 . Our results are in agreement with Zhang's argument, stating that 
the negative tail exponent rj(d) of the asymptotic behavior lnP(E'o) ~ — \Eo\ v ^ as Eq — > —00 is directly related to 
the fluctuation exponent 8(d) via the simple formula rj(d) = 1/(1 — 0(d)). 

Along the paper, we have also discussed the similarities and differences with spin-glasses. In particular, we have 
argued that the application of Zhang's argument for the Sherrington-Kirpatrick model of spin-glasses points towards 
an asymptotic distribution which is not a generalized Gumbel distribution g m {x), in contrast with the current way 
of fitting the numerical data Q, Hol | , but involves instead some non-trivial negative tail exponent tjsk > 1 directly 
related to the fluctuation exponent (Eq. I20|) . The fact that the fitting value m of generalized Gumbel distribution 
g m (x) depends on the probed range in the variable x (m ~ 6 via simple sampling |10| and m ~ 11 via importance 
sampling 0) also points towards t]sk > E More generally, we have explained in details how fits with generalized 
Gumbel distributions of the core of the distribution could be very misleading if one is interested on the tails, since all 
Gumbel distributions correspond to the exponent 77 = 1, which is very restrictive. 

Finally, our conclusion concerning the algorithm is that the importance sampling Monte-Carlo Markov chain in 
the disorder introduced in pj is a very efficient method to probe precisely the tails of probability distributions over 
the samples. In the field of disordered systems, this Monte Carlo procedure will be very useful to study probability 
distributions of other observables, beside the ground state energy. 
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